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ABSTRACT 


Previous work on the evolution of weak double layers in a hydrogen plasma has been extended to include H + 
and O with relative drift. It has been shown (Bergmann and Lotko, 1986) that the relative drift between hydrogen 
and oxygen ions due to a quasi-static parallel electric field gives rise to a strong linear fluid instability which 
dominates the ion-acoustic mode at the bottom of the auroral acceleration region. This ion-ion instability can 
modify ion distributions at lower altitudes and the subsequent nonlinear evolution of weak double layers at higher 
altitudes in the ion-acoustic regime. We have found that ion hole formation can occur for smaller relative electron- 
ion drifts than seen in previous simulations, due to the hydrogen-oxygen two-stream instability. This results in local 
modification of the ion distributions in phase space, and a partial filling of the valley between the hydrogen and 
oxygen peaks, which would be expected at higher altitudes on auroral field lines. It is shown that the observed 
velocity diffusion does not necessarily preclude ion hole and double layer formation in hydrogen in the ion-acoustic 
regime. These simulation results are consistent with the experimentally measured persistence of separate hydrogen 
and oxygen peaks, and the observation of weak double layers above an altitude of 3000 km on auroral field lines. 


I. INTRODUCTION 


Weak double layers with potential jumps comparable to the electron thermal energy have been observed to 
form in one-dimensional (Sato and Okuda, 1980) and two-dimensional (Barnes et al., 1985) electrostatic particle 
simulations; the double layer formation is driven by an electron drift relative to ions which is unstable to the ion- 
acoustic mode but is less than the electron thermal speed. Such weak double layers have been observed in space in 
the auroral particle acceleration region (Temerin et al., 1982), and in laboratory plasmas (Chan et al., 1984; Sekar 
and Saxena, 1985; Chan, 1986). Thus far, theoretical efforts at understanding weak double layer formation have 
focussed on a single ion species, while it is known from space observations that weak double layers occur in regions 
of upward flowing hydrogen and oxygen of ionospheric origin. A quasi-static parallel electric field has been postu- 
lated to explain the observed particle distributions (Chiu and Schulz, 1978; Lyons, 1980). While the existence of 
such a field will remain a zeroeth order assumption in the present paper, we will also examine non-adiabatic mod- 
ifications of the particle distributions at the bottom of the acceleration region which may affect weak double layer 
evolution further up the field line, and the stability of the assumed quasi-static field. 

We will first briefly review previous work on weak double layer evolution in a hydrogen plasma, and then 
extend our simulations to include a relative drift between hydrogen, oxygen, and electrons which occurs, for 
example, in a mirror-supported parallel electric field. Our purpose is to examine the nonlinear effects of the result- 
ing hydrogen-oxygen two-stream instability (Bergmann and Lotko, 1986) on the particle distributions, and con- 
sequences for double layer formation further up the field line. 


-'RECEDING PAGE BLANK NOT FILMED 


225 


III. WEAK DOUBLE LAYER FORMATION IN A HYDROGEN PLASMA 


Barnes et al. (1985) showed in a series of one- and two-dimensional, bounded and periodic particle simula- 
tions that weak double layers with potential jumps comparable to the electron thermal energy form when the system 
is driven by an electron drift relative to ions which is less than the electron thermal speed, e.g. , V d = 0.5 - 0.7 a e , 
and an electron to ion temperature ratio T e /T, >> 1 . The electron drift was maintained by injection of electrons from 
the boundaries at a continuous rate in bounded runs, and by applying a weak electric field uniformly across the 
system in periodic runs. 

Sato and Okuda (1980) first studied the occurrence of weak double layers in a one-dimensional periodic 
system in which electrons are given an initial drift that subsequently decays. They found that it was necessary to use 
a long system, L > 256 \ D (Debye lengths), in order for weak double layers to form in periodic runs. Our subse- 
quent interpretation (Barnes et al., 1985) is that long periodic systems are required to prevent electron recycling 
from the low to high potential side, which neutralizes the double layer. Electron injection boundary conditions 
eliminate this problem in bounded simulation runs, and a weak applied electric field acts to impede electron 
recycling in periodic runs; both of these simulation techniques allow shorter system lengths. 

Figure 1 from Barnes et al. (1985) shows the temporal evolution and recurrence of weak double layers in a 
one-dimensional system with electron injection boundaries. Ion-acoustic turbulence evolves, for V H = 0.5 a e and 
T e /Tj = 50, into a discrete localized pulse which propagates into the system initially at the sound speed. The pulse is 
characterized by a negative potential dip which amplifies by momentum exchange with reflected electrons (Lotko, 
1983; Chanteur et al., 1983); the asymmetic reflection of electrons results in a potential jump downstream. As the 
negative potential dip grows, it traps ions, slowing down the pulse via mass loading until an effective Bohm 
criterion for existence of the double layer potential jump is no longer met. The latter requires that ions flow into the 
high potential side at or near the sound speed (Chen, 1974), achieved here by motion of the pulse in the ion frame. 
The potential jump then decays and ion holes (Chan, 1986) or ion-acoustic solitons (Sato and Okuda, 1981) propa- 
gate away from the high potential side to seed new double layer formation. The decaying ion hole, still apparent in 
phase space, recoils backward as it moves downward through the ion distribution. 

Barnes et al . ( 1 985) examined the persistence of weak double layers in two-dimensional magnetized simula- 
tions. Electron injection boundary conditions produce one-dimensional double layers which are roughly uniform 
across the system in the direction perpendicular to B. To examine the transverse scale, a doubly periodic system 
with a weak electric field imposed uniformly along B was employed. The magnitude of the electric field was such 
that the corresponding potential drop across the system was less than the electron thermal energy, or eEo/T e = 
0.6/1 60 \ D . Figure 2 shows transverse localization of weak double layers for strongly magnetized electrons (o) ce /a) pe 
= 3 is the ratio of electron gyro to plas ma frequency ). The transverse dimension appears to decrease with increasing 
magnetic field strength, scaling with V\ D 2 + p s 2 , where \ D is the Debye length and p s is the ion gyroradius at the 
electron temperature. The parallel scale length remains the order of tens of Debye lengths, as in one-dimensionality. 
Ion-acoustic turbulence becomes homogeneous and does not evolve into localized weak double layers in weakly 
magnetized (a) ce /Wpe < 1) periodic systems. One therefore might expect to see such structures in the auroral acceler- 
ation region, but not, for example, in the solar wind. 
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IV. ION HOLES IN MULTIPLE ION SPECIES PLASMAS 


To the double layer evolution problem we would now like to add the effects of multiple ion species, H + and 
0 + , with relative drift. This introduces an important complication noted by Bergmann and Lotko (1986). A quasi- 
static parallel electric field produces a relative drift between ionosp heric hyd rogen and oxygen ions which have been 
accelerated through the same potential drop, such that V H /V 0 = VMq/M h = 4. This situation is fluid unstable for 
parallel propagating modes when the relative H + -O + drift exceeds a minimum, determined primarily by ion Landau 
damping, up to a maximum value that is less than about twice the hydrogen sound speed C s = VTe/M H . This 
indicates that the ion two-stream instability (for parallel propagating waves) will be confined to the bottom of the 
acceleration region, since at higher altitudes the relative drift will exceed the upper bound for instability. It is likely, 
although it has not yet been demonstrated, that obliquely propagating modes may still be unstable for drifts exceed- 
ing this upper bound. The growth rate for the ion two-stream instability is larger than that for typical (electron-ion) 
current-driven instabilities, and one might expect significant modifications of the hydrogen and oxygen dis- 
tributions to occur. In particular, the unstable ion two-stream waves have phase velocities lying between the 
hydrogen and oxygen distributions, and one might expect some quasi-linear filling, that is to say, formation of tails 
on the high and low velocity sides of oxygen and hydrogen, respectively. This quasi-linear filling could, in turn, 
affect the ion-acoustic instability and double layer evolution at higher altitudes, when ion drifts relative to electrons 
become a significant fraction of the electron thermal speed, as required for double layer formation in hydrogen 
plasma simulations. The instability analysis and simulations require knowledge or assumptions about the electron 
distribution in the region of interest. Bergmann and Lotko (1986) have integrated the electron distribution func- 
tions, F(vj|, v x ), in the Chiu-Schulz ( 1 978) equilibrium model of a mirror-supported electric field to obtain an effec- 
tive one-dimensional distribution, f(V||). These electron populations include precipitating magnetospheric 
electrons, primary and secondary backscattered electrons, and those electrons which are trapped between the mag- 
netic mirror below and retarding electrostatic potential above. At an altitude relevant to the ion two-stream 
instability, the bulk of ionospheric electrons has been retarded at lower altitudes by the potential drop which 
produces the relative ion drifts. The sum of the remaining electron populations, shown in Figure 3, is essentially a 
stationary Maxwellian with a precipitating electron tail. Also shown in the figure is a Maxwellian fit for the first 
three moments as described by Bergmann and Lotko. 

We would like to examine the spatial evolution of the ion distribution functions along the geomagnetic field 
line at various distances above the bottom of the acceleration region (nominally at an altitude of 2000 km in Chiu and 
Schulz, 1 978), including the ion two-stream unstable regime near the bottom on up to altitudes where the ion drifts 
become comparable to the electron thermal speed, and where double layers have been observed (>3000 km alti- 
tude). Since our computer resources limit the simulation system to lengths less than or the order of 1000 Debye 
lengths, we examine instead a temporal evolution problem which differs from the spatial evolution problem in at 
least one respect. In the spatial evolution case, the ratio of the H + /0 + bulk drift velocity is VMo/M h , as oxygen and 
hydrogen are accelerated to the same energy as a function of potential at a given altitude. In steady state at a fixed 
altitude there will be a continuous flow of oxygen and hydrogen whose drifts differ by a factor of 1 to 4, respec- 
tively, but the hydrogen and oxygen ions passing that altitude at a fixed time will not leave the bottom of the 
acceleration region simultaneously, since hydrogen flows up the field line faster. This follows from the relation 


e<{> = 1/2 M H V H 2 = 1/2M 0 V 0 2 


( 1 ) 


which holds at any given altitude where the potential is e<j>. Alternatively, in a simulation system evolving in time 
with a uniformly applied E„, the ion velocity varies as 
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which results in an H + -0 + velocity of M 0 /M H , rather than (Mo/M H ) l/2 . Furthermore, depending on the strength of 
the applied electric field E 0 , the ions may accelerate so rapidly that the upper limit on the relative drift for the ion 
two-stream instability may be exceeded before nonlinear saturation can occur. In such a case, we would not see the 
full effects of wave-particle interactions on the ion distributions. 

With these caveats in mind, we performed a series of one-dimensional electrostatic simulations using the 
particle code ES 1 (Birdsall and Langdon, 1984), in a periodic system of length 240 X D , using 1 6,000 hydrogen and 
1 6,000 oxygen ions and 32,000 electrons. We varied the uniform applied electric field from eE/Te = 0, 1 .2/240 \ D 
to 2.4/240 X. D and applied it only to the ions in order to simulate the approximately stationary electron Maxwellian 
(Fig. 3) through which the outflowing ions accelerate. We did initial value runs with V H = V D = O at t = O and 
runs which were initiated with V H and V 0 in the range where the ion two-stream growth rate peaks. 

Figure 4 shows the nonlinear evolution of the ion two-stream instability for initial drifts V H = 1 -2 C s and V 0 
= 0.3 C s and a uniform applied electric field eE/Te = 2.4/240 X D . The electron-to-ion temperature ratio is T/Ti = 
20 and the mass ratios are M H /M e = 50 and M q/M h = 8. The choice of drifts V H /V 0 = 4 is intermediate between 
the spatial evolution case where V H /V 0 = VM/Mh = 2V2 and the temporal evolution case where V H /V 0 = 
M 0 /M h = 8 for our mass ratio. Variations about this set of parameters are discussed below. One observes the 
formation of a localized fluctuation in the potential similar to that seen in the previously described (single ion) 
simulations at a time when the hydrogen drift relative to electrons is Vh = 0.2-0. 3 a e . This drift is smaller by a factor 
of 2 than in the single ion runs previously shown. The localized wave is a result of the nonlinear evolution of the ion 
two-stream instability which occurs at lower relative drifts (V h -Vq) than the current-driven, ion-acoustic 
instability. The potential pulse is subsonic in the ion frame, and appears to propagate with the ions out the right-hand 
boundary and re-enter on the left. Periodicity of the system allows one to see that the pulse is continuous from the 
right through the left boundary of an adjacent frame, since the pulse has not moved much from frame to frame. (The 
frames are separated in time by 60 w/ 1 .) A significant localized potential jump e<f>/T e > 1 develops, but does not 
persist as far downstream as in cases where the relative electron-ion drift is larger (Fig. 11). We therefore hesitate to 
call this structure a double layer when the system is in the ion two-stream unstable regime, although its features are 
very similar to those shown in Figure 1 , when translated to a frame in which electrons are stationary and ions drift. 
One sees trapping of hydrogen and oxygen on the sides of the distribution functions corresponding to the phase 
velocities of the (ion two-stream) unstable waves, namely the low velocity side of hydrogen and the high velocity 
side of oxygen. It seems appropriate to call this structure an ion hole. 

We observed ion hole formation in the ion two-stream unstable regime for a range of parameters summarized 
in Table 1 . The ion two-stream instability was observed over a broader range of parameters (Bergmann and Lotko, 
1986) than was ion hole formation, which apparently requires large amplitude waves and occurs only for sufficient- 
ly rapid linear instability. Recall that the ion two-stream instability is limited in duration as the electric field 
accelerates ions into and out of the range of linearly unstable drifts. Ion hole formation did not occur in runs 2-4 until 
the hydrogen bulk was accelerated to 0.2-0. 3 a c . In run 6, with no applied electric field but the same initial drifts as 
run 3, ion hole formation was not observed. In run 7, also with no electric field, but with initial drifts in the range 
produced by the electric field in run 3 at the time ion hole formation was observed, an ion hole forms. Sato and 
Okuda (1980) saw weak double layer formation in a system 256 X. D long but not in one 128 \ D long. Our system 
length of 240 \ D is marginally long enough to allow a double layer to form in the absence of an applied electric field 
before periodic electron cycling neutralizes the evolving double layer space charge. We also did a run (8) using a 
bounded one-dimensional electrostatic code, FDW1 (Lawson, 1984), with constant particle injection maintained 
by an external circuit and floating potential at both ends of the system, but with parameters otherwise the same as in 
run 7. Ion hole formation in runs 7 and 8 is comparable, as shown in Figure 5. 
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In order to address the temporal evolution question, we performed two runs (9 and 10) with two different 
values of the applied electric field, eEo/T e = 1 .2, 2.4/240 X D , no initial ion drifts, and periodic boundary con- 
ditions. Ion hole formation was evident but weaker for the larger electric field (run 9) than in the initial drift case 
(e.g. , run 3), and absent for the weaker electric field (run 10) when compared at time such that C s < V^-Vq < 2 C s . 
Some ion heating occurs in the initial value runs (9 and 10) before the relative drifts are comparable to the initial drift 
runs, i.e., optimum for ion two-stream instability. In initial value runs the hydrogen and oxygen ion distributions 
separate more quickly than in the case of spatial evolution, and so spend less time in the range of unstable relative 
drifts, V H - V 0 2 C s . The maximum growth rate of the two ion-stream instability is the order of ■y/<j 0 p L , — 10' 2 
(Bergmann and Lotko, 1 986, Fig. 3) for the k modes in our simulation system of length 240 \ D and grid size 0.5 \ D . 
Both initial value (run 9) and initial drift (run 3) cases remain in the range of unstable drifts Vh-V g ~ 1 .2-2 C s a 
number of e-folding times, but there appears to be some difference between initializing the system in the linearly 
unstable regime and evolving through it. In the auroral problem, one expects spatial evolution and weaker electric 
fields, discussed below, to separate the drifts more slowly relative to the linear growth time. 


The question arises whether ion heating by the ion two-stream instability, evident in Figure 4, will affect 
double layer evolution at higher altitudes where the ion drift relative to electrons is larger. Figure 6 shows the initial 
ion and electron distributions for runs 9 and 10. Figure 7 shows the same distributions for run 9 at the time an ion 
hole is beginning to form, while Figure 8 shows the same distributions at a later time when the hydrogen drift is 
becoming significant relative to electrons. Figure 9 shows the particle distributions in the weaker electric field case 
at a time when the drift is the same as Figure 7. We would conclude from this set of figures tht there is no major 
modification of the hydrogen and oxygen distributions by the ion two-stream instability, which is present in runs 9 
and 10. There is some heating on the low velocity side of hydrogen and the high velocity side of oxygen, as 
expected, in the range of unstable ion two-stream phase velocities. Figure 10, a similar plot for run 3 which showed 
ion hole formation at large trapping amplitudes (e<}>/T e ~ 1), indicates more heating between the hydrogen and 
oxygen distributions. This plot exhibits distribution functions which are spatially averaged across the whole 
system, and it is the case that the plateau evident in hydrogen (and oxygen) is due primarily to the spatially localized 
ion hole evident in Figure 3. It is questionable to call this heating versus localized ion trapping since it is not uniform 
across the system. It therefore seems reasonable to conclude that in our simulation system the hydrogen and oxygen 
average distributions are not so greatly modified by the ion two-stream instability as to preclude ion-acoustic 
instability and double layer formation at larger ion drifts relative to electrons. 


IV. ION-ACOUSTIC DOUBLE LAYERS IN AN 0 + -H + PLASMA 


As hydrogen and oxygen continue to accelerate up the geomagnetic field line out of the region of ion two- 
stream instability, hydrogen eventually acquires a drift relative to electrons comparable to the electron thermal 
speed. If the hydrogen velocity distribution has not been too greatly modified by the ion two-stream instability, as 
our preceding results indicate, we might expect double layers to evolve from the electron-hydrogen acoustic 
instability, as described by Barnes et al. (1985), as long as oxygen and hydrogen remain well separted. To test this, 
we did a series of simulations at large electron drifts (0.7-0.9ae) in the oxygen frame with hydrogen drifting at -0.2 
C s . Electron injection boundary conditions were employed as in Barnes et al. (1985), in a system 512 X D long, 
containing 8000 ions of each species and 16,000 electrons. More realistic mass ratios, M H /M e = 1000 and Mo/M H 
— 10, were used. Oxygen was kept cold relative to electrons, T e /T q = 100, and two cases were examined for 
hydrogen: T e /T H = 20 corresponding to no significant heating of hydrogen by the ion two-stream instability at 
lower altitudes, and T e /T H = 2 where significant heating has occurred. The assumption that electrons are hotter than 
ions is justified by the altitude where weak double layers have been observed (>3000 km), since the large scale 
parallel electric field restricts colder electrons to lower altitudes. No electric field was applied in these bounded 
runs. 
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Figure 1 1 shows the hydrogen and oxygen distributions and potential at a time when one and possibly a 
second double layer are forming with hole(s) evident in hydrogen phase space. Oxygen responds more slowly and 
appears to play a passive role in the double layer formation, but eventually forms a hole in ion-phase space by the 
time hydrogen has undergone significant heating and the double layer is disappearing (Fig. 1 2). In a similar run with 
T e /T H = 2, a hole does not appear to form in hydrogen but is evident in oxygen at later times. This result is con- 
sistent with Schamel’s( 1982) criterion that ion holes do not form forT e /Tj < 3.5 (see also Hudson et al., 1983). An 
oxygen ion hole and weak double layer appear to form when the hydrogen is heated too much to support such a 
structure. Should the hydrogen be significantly heated and the oxygen remain cool, a hole can still form in oxygen in 
association with an electron -0 + drift instability at phase velocities between the electron and 0 + peaks. 


V. APPLICATION TO THE AURORAL REGION 


A number of caveats are in order before applying the foregoing simulation results to the auroral particle 
acceleration region. We have examined separately two regimes: (1) where the two ion-stream instability operates at 
low relative ion drifts (<2 C s ) produced by a quasi-static parallel electric field near the bottom of the acceleration 
region, and (2) ion-acoustic double layer formation at higher altitudes where relative electron-ion drifts are larger 
and the ion distributions will have undergone some heating at lower altitudes. We have restricted our analysis to 
parallel propagating modes and one-dimensional simulations in the present paper. It is likely that oblique modes 
will affect the ion distributions. Kaufmann et al . ( 1 986) have examined the stability of hydrogen and oxygen beams 
measured by the S3-3 and DE 1 satellites and concluded that oblique modes are unstable. It may also be the case that 
the upper limit on relative drift for the ion two-stream instability is relaxed for oblique modes, since k • (V H - V Q ) <2 
C s for larger relative ion drifts when k is oblique. The linear instability of oblique modes is under investigation by 
Bergmann (private communication, 1986). Barnes et al. (1985) showed that ion-acoustic double layers evolve in the 
presence of oblique (EIC) modes in two-dimensionality with behavior similar to the one-dimensional case. We plan 
to extend the present multi-ion studies to two-dimensional in the future. 

Another qualification to our conclusions is the strength of the electric field used in the initial value simula- 
tions. The values of eEo/T e = 1 .2-2.4/240 correspond to 5-10 mV/m for == 10m, T e = 10eV,andn e = 10 
cm' 3 . These are not large parallel electric fields compared with observations in the acceleration region (Temerin, 
private communication, 1986), but are larger than the mirror-supported fields calculated by Chiu and Schulz ( 1 978) 
which maximize at E„ < 0.5 mV/m near the bottom of the acceleration region. 

It is somewhat difficult to extrapolate from the temporal evolution approach taken in this paper to the spatial 
evolution of ion distributions along auroral field lines. Nonetheless, with the neglect of oblique modes and use of 
somewhat large electric fields, and/or initializing the simulations with unstable drifts, and preheated ions in the ion 
acoustic regime, we find the following: 

1 . Ion holes form in the ion two-stream unstable regime at relatively low drifts compared with those requi- 
red to form ion-acoustic double layers. They occur in systems with and without an applied electric field, but over a 
narrow range of relative hydrogen-oxygen drifts. 

2. The ion two-stream instability does not appear to greatly modify the ion distributions, except locally in 
the presence of a large amplitude (e<t>/T e ~ 1) ion hole. 

3 . Double layer evolution should proceed at higher altitudes as previously described in the ion-acoustic 
regime, with noies forming in hydrogen, and oxygen responding passively. 
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There are two pieces of experimental evidence supporting our conclusion that hydrogen and oxygen dis- 
tributions are not so greatly modified at lower altitudes by the ion two-stream instability as to preclude double layer 
formation at higher altitudes.The first is the set of particle measurements from the S3-3 and DE 1 satellites analyzed 
by Kaufmann et al. (1986) showing well separated H + andO + peaks, with a slight filling in between the two. The 
second is the observation of what have been identified as weak double layers by Temerin et al. (1982), also Temerin 
and Mozer (1986), at altitudes >3000 km in regions of upward ion flows. These observations indicate that the 
heating of the ion distributions that occurs at lower altitudes is not as great as the relative acceleration, nor enough to 
make T e /Tj ~ 1 , which would preclude ion-acoustic double layer formation. 

A number of questions remain to be addressed: Can one design a temporal evolution simulation which better 
models the spatial evolution problem within the constraints of c omputer time and memory, e.g. , by moving one ion 
species with an electric field which differs from the other by VM H /M 0 to mimic the spatial case in an initial value 
problem? What effects do oblique modes and two-dimensionality introduce in the problem? Is it possible to use 
weaker electric fields and follow the evolution from the ion two-stream through the ion-acoustic double layer 
regime. These and other questions remain to be addressed in future work on the evolution of weak double layers in 
the multi-species auroral plasma. 

TABLE 1. ION HOLE FORMATION IN THE ION TWO-STREAM 
UNSTABLE REGIME 


Run 

V o (0) 

(3e) 

V H (0) 

(3e) 

Vh-V 0 

(C s ) 

eEo/Te 

be 

t 

(<v _ ') 

hole 

1 

0 

0.17 

1.2 

2.4 

per 

0 

no 

.2 

0.035 

0.17 

0.95 

2.4 

per 

0 

yes 

3 

0.042 

0.17 

0.91 

2.4 

per 

0 

yes 

4 

0.049 

0.17 

0.86 

2.4 

per 

0 

yes 

5 

0.06 

0.17 

0.78 

2.4 

per 

0 

no 

6 

0.042 

0.17 

0.91 

0 

per 

0 

no 

7 

0.06 

0.30 

1.7 

0 

per 

0 

yes 

8 

0.06 

0.30 

1.7 

0 

bnd 

0 

yes 

9 

0 

0 

0 

2.4 

per 

0 

marginal 

10 

0 

0 

0 

1.2 

per 

0 

no 


Note: Units of eE (> /T c are (240 X. D )‘‘; bnd and per refer to bounded and periodic boundary conditions (be). 
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Figure 1. Time series plots of ion-phase space (left) and electrostatic potential (right) for a bounded one- 
dimensional run with M/m = 2000, T/Tj = 50, and V H = 0.5 ae- The snapshots are taken at intervals of 360 tope' 1 
(8 oipj' 1 ) beginning at 1080 aipe' 1 (24 oipj 1 ) (from Barnes et al., 1985). 
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Figure 2. Magnetized doubly periodic system with (o ce /a) pe = 3 and a weak uniform applied electric field pointing 
to the left, eCo/Tg = 0.6/160 X d . (a) Potential profiles are averaged over 4 (Upf 1 ; (b) ion-phase space v x versus x is 
displayed in each of four bands in y. v x is in units of a e . The most prominent double layers are in band I at x ~ 1 20 \ D 
and band 3 at x =» 40 X D . 
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Figure 3 . Reduced electron distribution function f(V) for Chiu and Schulz ( 1 978) Model W, at e«J>/T c =110, where 
T c is the cold ionospheric electron temperature. This represents a sum of contributions from magnetospheric, 
primary and secondary backscattered and trapped electrons, with ionospheric electrons negligible at this value of 
the mirror potential measured from zero at an altitude of 2000 km (adapted from Bergmann and Lotko, 1986). 
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Figure 4. (Concluded) 











Figure 5 . (a) Potential for run 7 at w,*. At = 900; (b) Potential for run 8 at < 0 ^ At = 900. Parameters in (a) and (b) are 
the same, except that (a) was periodic, using ESI , and (b) was bounded using PDW1 , with particle injection main- 
tained constant as initialized (V H = 0.30, V Q = 0.06, V e = 0 in units of a e ) by an external circuit. All potential 
plots shown are averaged over 30 top,,' 1 . 
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